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Given n linearly independent pure states and their prior probabilities, we study the optimum 
unambiguous state discrimination problem. We derive the conditions for the optimum measurement 
strategy to achieve the maximum average success probability, and establish two sets of new equations 
that must be satisfied by the optimum solution in different situations. We also provide the detailed 
steps to find the optimum measurement strategy. The method and results we obtain are given a 
geometrical illustration with a numerical example. Furthermore, using these new equations, we 

■ derive a formula which shows a clear analytical relation between the optimum solution and the 
\ n states to be discriminated. We also solve a generalized equal-probability measurement problem 

^ . analytically. Finally, as another application of our result, the unambiguous discrimination problem 

O ■ of three pure states is studied in detail and analytical solutions are obtained for some interesting 

I cases. 

O \ PACS numbers: 03.67.-a, 03.65.Wj, 03.65.Ta, 42.50.Dv 
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O ^- I. INTRODUCTION 

^ ■ 

■ Discrimination of quantum states has been an interesting and attractive problem in quantum information science for 
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a relatively long time [1, 2, 3]. Since deterministic and error-free discrimination of an arbitrary set of quantum states 
is generally impossible due to the basic principles of quantum mechanics, just like other "no-go" theorems [1, H, 0|, 
the problem of finding an eflFective scheme of state discrimination has attracted a lot of attention and has played an 
import role in the study of quantum communication and cryptography. Considerable work has been devoted to this 
\ problem and it has developed rapidly recently. 

The task of state discrimination is to discriminate the state of a quantum system from a given finite set of possible 
states with certain prior probabilities, and there are mainly two kinds of strategies to complete this task. One kind of 
. strategy is called minimum error discrimination 0, 0, S| which requires that the average probabilityof identifying a 
■ wrong state is minimized. There have been numerous results for this kind of discrimination strateg y [3, IiqI. [ill 
H, [iH, , some of which are quite interesting, like the weighted square-root measurement flTl. fisl . flil [jol. [2ll. [22 
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■ and that measurement sometimes does not aid in discriminating certain set of states [25[, [26 
I Another important kind of discrimination strategy is unambiguous discrimination pioneered by [27|, |28|, |29| , which 
requires that no error occurs in the identification of the states at the expense of obtaining an inclusive result with some 
non-zero probability. A lot of research has also been performed on this kind of discrimination strategy [sO, [3l|, [H, [H, 
0, [H, [Sy , [13, [3^ , including unambiguous discrimination of symmetric states [39] and unambiguous discrimination 
between mixed states [10, [i3, [H, [13, H, [H, [H, [13, [H, [li] . 
\ In addition to the above two strategies, research has also been performed on mixed strategies involving minimum- 
error discrimination and unambiguous discrimination together [HO, [HI, [H, [H^l , in order to achieve balance between the 
accuracy and the eflftciency of state identification. It is, in general, very diflftcult to maximize the average probability of 
successfully discriminating the given states analytically, but some special techniques such as semidefinite programming 
have been employed to solve this computation problem numerically (ssl. IH^. IHBI. [56| . 

In this article we shall study the properties of the optimum strategy for unambiguous discrimination of n pure 
states, and give a detailed method to obtain such an optimum strategy. According to (31I . [32| . one can manage 
to discriminate a set of states unambiguously with non-zero success probabilities if and only if the given states are 
linearly independent. We shall hold this assumption throughout this article. The main method we use for studying 
this problem is to put the individual probabilities that each state is successfully identified together as a vector in the 
n-dimensional real space (and do the same to the n prior probabilities), and study the properties of the minimum 
eigenvalue of the matrix X — V under the optimum strategy by vector analysis techniques (X is the Gram matrix of 
the n states to be discriminated and V is the diagonal matrix with the success probabilities as its diagonal elements, 
they will be defined explicitly later in Theorem [1]). From the properties of the minimum eigenvalue of X — F, we shall 
establish two sets of equations that the optimum solution of the unambiguous discrimination problem must satisfy in 
diflFerent situations, and these equations will turn out to have some intuitive geometrical meanings. 
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This article is organized as follows. In Sec. [Til we shall give a general description of the unambiguous state 
discrimination problem that we are interested in and the POVM formalism that we shall use in this article. In Sec. 
IIIH we derive the properties of the optimum measurement strategy and the new equations that can be used to work 
out the optimum solution. Examples are also given to geometrically and numerically illustrate our method of solving 
the unambiguous discrimination problem. Sec. IIVI is devoted to deriving an analytical formula which characterizes 
a simple relation between the maximum average success probability and the n pure states to be unambiguously 
discriminated. In Sec. [Vl a generalized version of equal-probability-measurement (EPM) problem [HHl is studied and 
an analytical solution is obtained. Finally we apply the results obtained in Sec. IIIII to the case of three linearly 
independent pure states in Sec. IVII and work out analytical solutions for some interesting cases. 



II. PROBLEM DESCRIPTION AND POVM FORMALISM 



In this article, our problem is how to unambiguously identify the state of a quantum system with the maximum 
average success probability since the discrimination of non-orthogonal states is generally probabilistic. And all we 
know is that this state belongs to a given set of n linearly independent states with given prior probabilities 

7i (z = 1, • • • , n). We would like to obtain some analytical conditions that the optimum solution should satisfy, and 
provide a detailed method for obtaining the optimum solution of the problem. 

In the following study of the problem, the n prior probabilities will be denoted by a real vector 7 in the space 
and the success probabilities (z = 1, • • • , n) for unambiguous outcomes will also be denoted as a real vector p in 
for short. The average success probability for unambiguous discrimination can be written as p = Xir=i 7^^^ —IP 
then. 

Since only the probabilities of the measurement outputs are concerned in this problem, we shall use the POVM 
{positive- operator-valued measure) formalism [57] which is a good description for the statistics of a general physics 
process. A POVM consists of a set of POVM elements {11%} satisfying 11^ > and 11^ = /, where / represents the 
identity operator (or matrix). A Hermitian operator is said to be positive and denoted by "> 0" if all of its eigenvalues 
are non-negative. To unambiguously identify the state, we require that 

P{j\i) = iH^jl^i) =PAj^ Vz,j = I,-- - ,n (1) 

where p{j\i) represents the probability of obtaining the result j when the original state of the system is actually 
and Pi denotes the probability of correctly identifying {tpi). 

Let H denote the n-dimensional Hilbert space spanned by the given set of linearly independent pure states 
The total Hilbert space of the system, denoted by Htotai = H ® 7Y^, may be larger than 7Y, where denotes the 
subspace orthogonal to Ti. Any operator acting on Htotai is equivalent to an operator on 7i by projecting onto the 
space H when only the effect on H is considered. Therefore, without loss of generality, we shall restrict our POVM 
elements to those acting on the space H spanned by the given set of states. In order to discriminate the states 
{\ipi)}^^i unambiguously, the POVM element Hi that identifies the ith state {ipi) must be orthogonal to the subspace 
spanned by the other n — 1 states according to Eq. ([T]) , therefore the rank of Hi should be no larger than 1, for all 
z = 1, • • • , n. Thus each POVM element that successfully identifies a certain state should have the form [sl, [HHl 

n, =p,|v^,)(V^,|, (2) 

where pi is the success probability to identify the ith state and \ipi) is an unnormalized state orthogonal to {t/jj) for 

all j ^ i. It can be seen by substituting Eq. ([2]) into ([T]) that [(^/^il^/^i) | =1 and there is freedom for each li/ji) to have 
an arbitrary phase, and without loss of generality we make a specific choice of the phase such that 

(ijjli^i) = Sij, (3) 

for convenience. _ 
Let ^ denote the matrix with Itpi) as its ith column. Define <l> as 

$ = <|>(<|>t<|>)-i, (4) 

then 

<l>^i = /. (5) 

Note that Eq. dU cannot be simplified to ^ = (^^)~^ in general, because the states may belong to a larger 

Hilbert space of which the dimension is greater than n, implying that the matrix <I> may not be a square matrix. (In 



3 



the last paragraph, we only restrict the POVM elements to be those that act on the space spanned by these states, 
but we do not make any restrictions on the representation of these states.) Comparing (|3]) and (O, we know that 
is exactly given by the ith column of the matrix <l>. 

Since unambiguous discrimination of the states {{ipi)}^^-^ is probabilistic if they are not orthogonal to each other, 
there exists a POVM element Ho which gives the inconclusive result, and it can be written as 

n 

^o = I-J2pM^^m (6) 
Since any POVM element must be positive to represent a physically realizable process, it is required that 

n 

no = /-^]?#.)(^?i| >o. (7) 

This positivity inequality is an essential constraint on the unambiguous discrimination scheme and is the starting 
point of the discussions in this article. 

When the set of states and prior probabilities 7^ = 1, • • • , n) are given, the POVM elements for the 

measurement {Hi = pi\ipi){tlji\, Hq = I — X^ILi depend only on the variables pi, i.e. the success probabilities, since 
l^pi) (i = 1, • . • , n) can be determined and explicitly given by the ith column of ^ defined in (|4]). Therefore, searching 
for the optimum solution is to find a set of success probabilities Pi (corresponding to a point or a vector p in R^) such 
that their weighted average (with the prior probabilities 7^ as the weights) is maximized under the restrictions Pi > 
(i = 1, • • • , n) and Uq > 0. The optimum solution is denoted by Popt^ shall be called optimum point sometimes 
throughout the rest of the article; and the main goal of this article is to find the optimum point Popt- 

III. EXPLORATION OF THE OPTIMUM STRATEGY 

In this section, we shall obtain some properties of the optimum strategy for the unambiguous state discrimina- 
tion problem, and provide a systematic way to obtain the maximum average success probability and the optimum 
measurement strategy. 

A. General properties and methods 

As the first step to study the optimum unambiguous state discrimination problem described in Sec. [TTl we are 
going to re-derive the positivity condition given by Duan and Guo [32| in a more concise way, using the POVM 
representation, and prove a convexity property of the set of all feasible p's. 

Theorem l(Positivity and Convexity). Suppose {IV^i)}^^^ ^ of linearly independent pure states. Let X = ^"^^ 
where ^ is the matrix whose ith column is (i = l,---n), and F = diag(pi,--- ,Pn) where pi is the success 
probability to unambiguously discriminate 1?/;^) (i = 1, • • • , n). Then i) 

X - F > 0, F > 0. (8) 

ii) Let S denote the set of points p satisfying the positivity condition (|8]), then S is convex. 
Proof, i) Using the definition of F in this theorem, Eq. ([7j) can be re-written as 

/ - iF^^ > 0. (9) 

Substituting Eq. (|4|) into (|9]) , we have 

/ - $(<i>t$)-ir($t<i>)-^^t > 0. (10) 

According to the property of positive matrix, we can multiply Eq. ([TQ|) by <l>^ from the left side and by <l> from right 
side, therefore immediately get the first inequality in Eq. (|8]). The second inequality of Eq. (|8]) must be satisfied since 
the success probabilities > (i = 1, • • • , n) must be non-negative. 

ii) Let Pi and P2 denote two arbitrary points in 5, and Fi, F2 denote the diagonal matrices with the components 
of Pi and P2 as their diagonal elements, respectively, then X — Fi > 0, Fi > 0, X — F2 > 0, F2 > 0. Let e be an 
arbitrary real number between and 1 and 



Pe = + (1 - e)p2, 



(11) 
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re = eri + (i-e)r2. (12) 

Since the sum of two positive matrices is still a positive matrix, we have 

X-T,=X- (ePi + (1 - e) T2) = e (X - Ti) + (1 - e) (X - T2) > 0, (13) 

re = eri + (i-e)r2>o, (i4) 

which means that e S. Thus 5 is a convex set. ■ 

Before deriving more properties of the optimum solution for unambiguous discrimination, we introduce some nota- 
tions and nomenclatures that will be used later on here. 

Let (Ti, • • • , denote the eigenvalues of X — F in decreasing order such that cri > • • • > cr^, then the minimum 
eigenvalue must satisfy cr^ > due to the positivity of X — F according to Theorem [TJ Since X — F depends on 
the parameters {pi, • • • ,Pn}, which is denoted by a real vector p for short, we have = cr^(p), (i = 1, • • • , n). 

For the convenience of description, the following notations and nomenclatures will be used throughout this article. 

• a point: a vector p in R^; 

• the feasible set: the set S defined in Theorem [U 

• the critical feasible region (denoted by IZcf)'- the set of points p in S satisfying an {p) = and F > 0; 

• the boundary of the critical feasible region (denoted by Br): the set of points p in the critical feasible region 
1ZcF with at least one zero component, i.e., Pi = for at least one i G {!,•'• ? 

• the interior part of the critical feasible region (denoted by ftR): the set of points p in the critical feasible region 
IZcF with Pi > for alH G {1, • ' ' 7 ^l? i-^. = ^cf\Br where "\" denotes the set exclusion operator; 

• an interior point: a point in the interior part VLr of the critical feasible region (but not off the critical feasible 
region) ; 

• a boundary point: a point on the boundary Br of the critical feasible region; 

• a singular point: a point in the critical feasible region where Vcr^ (p) does not exist (i.e. cr^ (p) is degenerate) 
or Vdn (p) = 0. Here V denotes the gradient operator. 

It is clear that Br U Q.r = 1Zcfj and 1Zcf C S. 

We give an example of three pure states below to explain the nomenclatures defined above intuitively by graphics. 
Suppose the states to be discriminated are 

1^1) = (1, 0, of , 1^2) = ^ (1, 2, of , 1^3) = ^ (1, 1, (15) 

Fig. [1] shows the feasible set 5, the critical feasible region 1Zcf, the boundary Br of the critical feasible region and 
the singular point pg for these three states. But it should be noted that the singular point does not always exist for 
every set of states. A necessary condition for the existence of a singular point in the case of three states is given in 
Case 3 of Sec. IVIBi In addition. Fig. [1] explicitly shows that the feasible set is convex, which agrees with Theorem 

m 

Now we give two lemmas that will be used in the proofs of the next several theorems. 

Lemma 1. If there exists a point pg ^^le critical feasible region IZcf satisfying that ^ ■ p < 'J ■ Po for any point p 
in a sufficiently small neighborhood A of Pq, A C S, then 7 • p < 7 • Po holds for any point p G <S. 

Proof. By contradiction. Suppose there exists such a point p^ G 5 that 7 • Pi > 7 • Pg, let e be a sufficiently small 
positive real number satisfying < e < 1, and let 

Pe = ePi + (1 - e)Po, (16) 
then p^ ^ S because of the convexity of S proved in Theorem [H and p^ G A since e is sufficiently small. However, 

7 • Pe = 7 • (ePi + (1 - e) Po) = ^7 • Pi + (1 - e) 7 • Po > 7 • Po (1'^) 
which contradicts the assumption of this lemma. Therefore, 7 • p < 7 • Po holds for any point p in SM 



5 




Figure 1: (Color online) Feasible Set, critical feasible region, boundary and singular point. The whole shaded volume (including 
its surface) is the feasible set <S, the upper curved surface of the feasible set is the critical feasible region IZcf^ and the black 

curved boundary line BC U AC U AB (AB is hidden and invisible in the figure) is the boundary Br of the critical feasible 
region. The singular point is the point at which the curved surface is not smooth. 



Intuitively, Lemma [T] tells us that if a linear function acting on a convex set achieves a local maximal value at some 
point of the set, then that local maximal point must be the global maximum point that the function can reach over 
the whole set. 

Lemma 2. Suppose a and b are two non-zero real vectors in the space, {a - x) {b - x) < for any vector x, if 
and only if a and b are anti-parallel, i.e., a = —Xb where A is positive factor. 

Proof. Necessity: By contradiction. Let x = ao + 6o, where «o = and 6o = j^y. Let denote the angular between 
a and 6, i.e., a - b = \a\\b\ cos6. If a and b are not anti-parallel, then < 6 < n hj the well known Cauchy- Schwartz 
inequality [5§| and x 0. Then we have (a ■ x) {b ■ x) = \a\\b\ (1 + cosO)'^ > 0, contradictory to the assumption. 

Sufficiency: If a and b are anti-parallel, i.e., a = —Xb (A > 0), then for any non-zero x, (a - x) {b ■ x) = — A (6 • x)^ < 

0. ■ 

Now we are going to show some important properties of the optimum point, from which we will establish two sets 
of new equations that the optimum solution has to satisfy later in different situations. 

Theorem 2 (Minimum Eigenvalue), i) The optimum point p^^^ must be a point in the critical feasible region TIcfj 

1. e., 

{Popt) = 0. (18) 

ii) If the optimum point p^p^ is a non-singular point in the interior part Qr of the critical feasible region, p^^^ G ftR, 
then 

V^n (p) |p,,, = -7. (19) 

iii) Conversely, if there exists a point satisfying Va^ (p) Ip^ = —7 in the critical feasible region TZcfj Po must be 
the global optimum point. 

Proof, i) Since cr^ is the minimum eigenvalue of the matrix X — T and X — T is positive as proved in Theorem [TJ 
we have an > 0. Thus, in order to prove Eq. ([18]), we only need to show that an cannot be greater than zero at the 
optimum point. We prove it by contradiction. 

If an ipopt) > 0, we let p^^^ be changed a little by then we can always find such 5p that 7 • (5p > and at the 
same time an {Popt + ^P) > 0, F + (5r > 0, where ST is a diagonal matrix with Spi^ • • • , Spn as its diagonal elements 
(such Sp always exists as we can construct it simply by taking Spi > and Spi sufficiently small, for alH = 1, • • • , n). 
This implies that the point p^^^ + Sp in the feasible set S satisfies 7 • {Popt + Sp) > 7 • Popt, which contradicts the 
assumption that 7 • p reaches the maximum at p^pt, so Eq. ([T8|) holds. 

It seems that the above proof may be applied to other eigenvalues of X — F in a similar way, but if any eigenvalue 
other than the minimal one is equal to zero, then the minimal eigenvalue would be negative and X — T would not 
keep positive, violating Eq. ([8]). 
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ii) Since p^p^ is an interior optimum point in the critical feasible region IZcf, if we change p^p^ a little by any 
Sp such that (Jn{Popt + ^P) ^ 0? must be that ^ - 5p < 0, otherwise 7 • p is not maximal at Pop^ Considering 
(Jn {Popt) — inequality crn{Popt + ^P) — ^ can be converted to 

5cJn ip) |p„,, > 0. (20) 

Since p^p^ is not a singular point, Van (p) \p^^^ exists, so 

Sdn (p) = (Jn {Popt + ^P) - {Popt) = (p) |p^^^ • Sp. (21) 

Substituting Eq. ([21|) into ([20]) , we have 

V^n (p) |p,^, • > 0. (22) 

Comparing Eq. ([22]) with 7 • (5p < 0, it follows that 7 and Van (p) Ip^^^ must be anti-parallel according to Lemma [2 
so 

Van (p) |p^^, = -^7, a > 0. (23) 

We will prove that a = 1 in the proof of Theorem [31 thus Eq. ([T9|) holds. 

iii) If there exists a point pg where Van (p) Ip^ = —7 in the critical feasible region 1Zcf, any small Sp such that 
(Po + ^P) = ^<^n (p) = Vdn (p) ' ^^P > will lead to 7 • (5p = —Van (p) Ip^ • ^^P < 0, so 7 • p reaches maximal 

at Po in a small neighborhood of pg, and according to Lemma [H Pq must be the global optimum point where 7 • p 
reaches the global maximum. ■ 

It should be pointed out that a result equivalent to our result i) in Theorem [2] was obtained in Ref. [31], where 
the maximum eigenvalue of <l>r <!>''' is proved to be exactly 1 when the unambiguous state discrimination scheme is 
optimum. Besides, Eq. ([T8|) implies that 

det (X - F) = (24) 

at the optimum point, and a physical interpretation of this fact was given in (33 |. 
To prove the next theorem, we first give another lemma [59]. 

Lemma 3. If ri,--- ,rn are the eigenvalues of an n x n matrix A, then ^ ^ii'''^ik — 

l<ii < • • •<ifc <n 

(all k X k principal minors of A). 

Based on Theorem [2 we now present the following theorem which can be conveniently used to obtain the maximum 
average success probability and the optimum measurement strategy when the optimum point is not singular and in 
the interior part Q.r of the critical feasible region. 

Theorem 3 (Interior Non-Singular Solution). Let Mk (p) denote the principle minor of order n — 1 associated with 
the kih. diagonal element of X — F, then a non-singular point pg in the interior part VLr of the critical feasible region 
is the optimum point if and only if pg is a solution of the following set of equations 

Mi{p) =7iA 

• , (25) 

Mn (P) = 7nA 

^det (X -F) = 

for some positive number A > 0, and satisfies the positivity condition ([8]) at the same time. 

Proof. First, we prove the "only if" part of the theorem. Let's consider the variation of det(X — F). Since 
det (X — F) = cTi • • • cTn, we have 

S det (X - F) = {Sai)a2 • • - an + ai {Sa2) ^3 • • - an H h aia2 • • -crn-i^an. (26) 

If Po is the optimum point that achieves the maximum average probability, then <Jn (Po) = according to Theorem 
[21 Eq. ([Ml) can be simplified to 

(5det(X-F)|p^ = (ai(72---an-i^an)|p,. (27) 

Since Pq is not a singular point. Van (p) Ip^ exists and San{p)\pQ = Van (p) Ip^ • Sp. From Jdet(X — F) = 
V det (X - F) • Sp and Eq. ([27|), we have 

V det (X - F) = aia2 ■ ■ • an-iVan (p) |p,. (28) 
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Again using (Jn (po) = 0, we can get 

V det (X - r) = (CTI • • • (Jn-l + CTi • • • an-2Crn H h Cr2 • • • CTn) Vdn (p) Ip^ • (29) 

According to Lemma [31 we have 

V det {X - r) = (Ml (po) + • • • + (po)) Va, (p) (30) 

On the other hand, 

Vdet (X - r) = f + • • • + ) det (X - T) , (31) 



dpi dpn 

where Ci is the orthonormal basis vector in the space associated with the coordinate pi. By performing Laplace 
expansion on the determinant of X — F along the kth row (or column) , we can see that 

^det(X-r) = -Mfc(p), (32) 

OPk 

so 

Vdet (X - F) = - (Ml (p) ei + • • . + Mn (p) e^) . (33) 
Comparing Eq. ([30]) and ([33]) and substituting Eq. ([23]) into ([30]), we have 

- (Ml (po) ei + • • • + Mn (po) en) = - (Mi (po) + • • • + M, (po)) ^7, (34) 

which results in 

Mk (Po) = o^T/c (Ml (po) + • • • + Mn (Po)) , /c = 1, • • • , n. (35) 

Considering that X — F > 0, we have (Po) > for all /c = 1, • • • , n. Since the theorem assumes that the optimum 
point is not singular, the zero eigenvalue is not degenerate and Van (p) {p^ 7^ 0, then according to Eq. ([28]) . we have 

Vdet(X-F)|p^ y^O. (36) 

So according to Eq. ([33]) we can see that the M^ (Po)'s are not all zeros and 

Ml (po) + • • • + Mn (po) > 0. (37) 

By summing up Eq. ([35]) for i = 1, • • • , n and using 71 + • • • 7^, = 1 and Eq. ([37|) . we immediately have a = 1. (This 
completes the remaining part of the proof for part ii) of Theorem [2) 
Thus, Eq. ([35]) can be simplified to 

Mk (Po) = Ik {Ml (po) + • • • + Mn (po)) , /c = 1, • • • , n, (38) 

which immediately leads to Eq. ([25]) . Considering Eq. ([37|) . there must be A > according to Eq. ([25]) . 
This concludes the proof of the "only if" part of the theorem. 

Next, we prove the "if" part of the theorem. Considering X — F > 0, det (X — F) = in Eq. ([25]) and that an (po) 
is the minimal eigenvalue of X — F, it follows straightforwardly that 

CTn (Po) = 0. (39) 

And since A > 0, Mk (Po) > 0, /c = 1, • • • , n, thus the rank of X — F is n — 1 and V det (X — F) 7^ 0, which imply 
that Po is not a singular point. It can be seen that Eq. ([30]) still holds here, so substituting ([25]) into Eqs. ([30]) and 
([33]) . one can have 

- (7iAei H h 7nAen) = AVan (p) Ip^, (40) 

where 71 + • • • 7n = 1 has been used. Therefore we get Va^ (p) Ip^ = —7, and together with ([39]) it implies that Po 
is the (global) optimum point according to part iii) of Theorem O This concludes the proof of the "if" part of the 
theorem. ■ 
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Theorem [2] and Theorem [3] describe the properties of the optimum point when it is a non-singular point in the 
interior part of the critical feasible region, and Theorem [3] also gives a way to find such an optimum point. When 
one obtains a solution from Eq. ([25|) for some A > 0, he or she has to verify whether it satisfies Eq. (|8]). If it does, this 
solution is the optimum solution. However, if ([25]) has no solution satisfying A > and ([8]), then the optimum point 
is either a point on the boundary Br or a singular point (if it exists) of the critical feasible region IZcf- The next 
theorem is to characterize the properties of the optimum point when it is on the boundary Br and give a method to 
work out the optimum point in that situation. 

For the simplicity of later description, we further define some new notations here. We denote the part of the 
critical feasible region IZcf where Pi-^ = 0, • • • , =0 as Br (ii, • • • , ife), and any Br (n, • • • , ik) will be called a 
{n — k — 1)- dimensional boundary of the critical feasible region (since the critical feasible region itself is of dimension 
n — 1 in the space). We again take the three state ([T5|) as an example. In Fig. [U the black boundary line excluding 
the points A, B and C is the 1-dimensional boundary of the critical feasible region, and the points A, B and C are 
the 0-dimensional boundaries of the critical feasible region. 

Theorem ^(Boundary Solution). A point pg on a (n — /cq — l)-dimensional boundary Br (n, • • • ,Zfeo) but not on 
any lower dimensional boundary, i.e., = 0, • • • , pi^^ = and pj > 0, V j 7^ n, • • • ^ikoj is the optimum point if and 
only if it is a solution of 

f^i(p)|p^i=o,-,p,,^=o = 7iA VzG {!,••• ,n}\{n,--- ,z/eo} 

|det(X-r)U^^o,..,p.,^=o = ' ^ ^ 

satisfying A > 0, the positivity constraints (|8]) and 



(42) 



Mi. (p) |p,^=o,-,p,, =0 > A7i, 



Proof. First, we prove the "only if" part of the theorem. If Pq is an optimum point on the {n — ko — l)-dimensional 
boundary Br (n, • ' ' ? ^/co) but not on any lower dimensional boundary, let Pq be changed to Pq + 6p on Br (n, • • • ^iko) 
by any small Sp satisfying Spi^ = 0, • • • , Spij^^ = 0, then there must be 

,5 det {X - r) |p„ = V det (X - r) • 5p = - (Po) ^Pi = 0, (43) 



where we have used Eq. (|33|) . and 



1 •Sp= ^ 7i5pi < 0. (44) 



So 

(Ml (po) , ■ ■ • , Mi (po) , • • ■ , M„ (po)) °^ (71, • • • (45) 

according to Lemma [21 therefore, Eq. ([4T]) holds. 

On the other hand, let pg be changed to pg + ^P' ^^e critical feasible region 1Zcf by another arbitrary small Sp' 
satisfying that Sp[^ > 0, • • • , > 0, then there must be 



,5det(X-r) \p^ = -Y,M,{po)5p[ 

1 

^ M,{p^)6p[- Yl M^iPo)^P^ (46) 



where we have used Eq. ([4T]) . and 



7-Sp' = Y^M= E E ^M<0- (47) 

2=1 i=iij---,ikQ 
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Substituting Eq. ([46]) into (|47l) , we have 

^.Sp' = j2^M= E U-^^)¥.<0. (48) 

Considering Sp[^ > 0, • • • , 6p[^ > and A > 0, it can been seen that Eq. (|42]) holds according to Eq. (|48|) . This 
concludes the proof of the "only if" part of the theorem. 

The "if" part of the theorem can be directly proved by reversing the above reasoning and using Lemma [H so we 
are not going to show the details here. ■ 

Since the feasible set 5 is a closed convex set confined in a finite region of R"^, the optimum point where the average 
success probability 7-p reaches the maximum always exists. Thus if Eq. (|4T]) for any {n — ko — l)-dimensional 
boundary (1 < ko < n — 1) does not have a solution satisfying A > 0, the positivity constraints (|8]) and Eq. (|42|) 
while the optimum point is also not an interior non-singular point in the critical feasible region IZcf, the optimum 
point can only be a singular point then. We know that a singular point is a point in the critical feasible region where 

(p) = is degenerate or Va^ (p) = 0, so when the optimum point is singular, Vdet {X — T) = according to 
Eq. ([28|) . which, together with (|33|) . implies that all (p) = for all /c = 1, • • • , n. Thus a singular point can be 
obtained as a solution of ([25]) with A = and the positivity conditions (|8]) . 

Remark 1. It should be pointed out that given a set of linearly independent states, if there exists a singular point 
in the critical feasible region IZcf, then that singular point could be the optimum point for a range of different 7's, 
since the normal vector of the critical feasible region IZcf changes discontinuously in the neighborhood of a singular 
point. 

In this subsection we have mainly studied the properties of the optimum point and obtained the equations that the 
optimum point must satisfy in different situations. We summarize our method to find the optimum point as follows: 

Step 1. Try to solve Eq. ([25]) in Theorem[3]and see whether there exists a solution satisfying the positivity constraints 
A > and Eq. ([8]). If such a solution exists, it is exactly the optimum point we try to find. 

Step II. If the set of equations ([25]) does not have a solution that satisfies A > and Eq. ([8]), one has to continue 
to search for the optimum point on the boundary Br of the critical feasible region using Theorem [H In detail, one 
can first solve Eq. (|4T]) on all (n — 2)-dimensional boundaries and see whether there exists a solution that satisfies 
A > 0, the positivity constraints ([8]) and Eq. (|42]) . If such a solution exists, it is exactly the optimum point; otherwise, 
one should further search on all (n — 3)-dimensional boundaries, (n — 4)-dimensional boundaries, . . . , until such a 
solution is found or all boundaries of dimension lower than n — 1 have been searched. If such a solution is found, it 
is exactly the optimum solution. 

Step III. If the optimum point is not found in the above two steps, then it must be a singular point and can be 
obtained by solving Eq. ([25]) with A = 0. If there exists more than one singular point on the critical feasible region 
IZcF, the one that maximizes 7-p is the optimum point. 

By Theorem [3] and [H we give two sets of explicit analytical equations for solving the unambiguous discrimination 
problem in different situations. We can use them to work out analytical solutions or obtain some analytical relations 
for the problem (see examples in Sec. \W\ (Vl and IVl]) . However, since Eqs. ([25]) and ([4T]) are nonlinear and the 
variables pi, • • • ,Pn are tightly coupled in the equations, maybe only numeric solutions can be obtained for these 
equations in some situations. A lot of sophisticated numerical techniques like Newton's method (including many of 
its variants), hybrid Krylov methods and so on have been developed to solve such nonlinear equations (gqI. IgiI. IgI. [63| . 

It is worth mentioning that other methods such as semidefinite programming [54 i55l i56i | have been developed to 
solve this unambiguous discrimination problem. Those methods are developed from some classical numerical analysis 
theories and they are mostly suitable for finding numeric solutions, while our method is developed purely by algebra 
and go in a totally different way, aiming at providing a new tool to treat the problem analytically. 



B. Geometrical view and numerical example 

A geometrical method is given in [s^l to solve the optimum unambiguous discrimination problem for three pure 
states, mainly for the situation when the optimum point is an interior non-singular point. In this subsection, we are 
going to give a similar but more complete geometrical way to illustrate the problem and the results we obtain in the 
previous subsection, and we calculate a numerical example illustrated with corresponding graphics to explicitly show 
the geometrical meanings. 

Geometrically, p = 7 • p = 71^1 + • • • + jnPn can be perceived as an (n — l)-dimensional plane in the space, and 
the critical feasible region an (p) = can be perceived as a curved surface in W^. 
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Table I: Numerical results of optimum points with different prior probabilities 



Category 




Prior Probabilities 




Pi 


P2 


P3 


A 


Popt 


Interior Point 


71 = 


= 0.05,72 


= 0.35,73 = 


0.60 


0.5029 


0.3169 


0.3629 


0.2326 


0.3538 


Boundary Point 


71 = 


= 0.10,72 


= 0.80,73 = 


0.10 


0.3927 


0.5300 





0.6577 


0.4632 


Singular Point 


71 = 


= 0.30,72 


= 0.35,73 = 


0.35 


0.6667 


0.4000 


0.2941 





0.4429 



It can be shown easily that the vertical distance from the origin of the coordinate system (with pi , • • • , Pn as the 
coordinates) to the plane p = 7 • p is 

(49) 



V7? + --- + 7^ 

so the average success probability p characterizes the vertical distance between the origin and the plane p = 7 • p in a 
geometrical view. Therefore, the problem of optimum unambiguous discrimination of pure states can be translated to 
the problem of finding an optimum point in the feasible set S at which the plane with fixed normal vector (71, 72, 73) 
(unnormalized) is most distant from the origin. Obviously the optimum point must lie in the critical feasible region 
TZcFj which is the "surface" of the feasible set 5, and this is in accordance with i) of Theorem [2l 

Now suppose the plane p = 7 • p is moved by parallel shifts, i.e., by changing p while the normal vector (71, 72, 73) 
keeps fixed. If the plane can be tangent with the critical feasible region IZcf when p is equal to some Pq, the distance 
from the origin to the plane is then maximized and the tangent point is exactly the optimum point PoptJ and po is 
the maximum average success probability. 

When the plane and the critical feasible region IZcf are tangent, their normal vectors at the tangent point should 
be parallel or anti-parallel. This implies that V (7 • p — p) = CV det (X — F), where C and C / 0- Using Eq. (|33|) . 
this equation can be simplified to Eq. ([25|) . 

However, if the plane can never be tangent with the critical feasible region when the plane is moved by any parallel 
shift, it means that a non-singular optimum point does not exist in the interior part of the critical feasible region, 
then the optimum point is either a boundary point on Br or a singular point in the critical feasible region. And if 
the optimum point is on some {n — k — l)-dimensional boundary, the plane is then tangent with that {n — k — 1)- 
dimensional boundary, resulting in Eq. (|4T]) . Conversely, if the plane is tangent with some {n — k — l)-dimensional 
boundary, it does not imply that the tangent point must be the optimum point though, unless Eq. (|42]) is satisfied, 
which ensures the average success probability at any other point in the feasible set will be no larger than that at the 
tangent point, due to the convexity of the feasible set. This gives the geometrical meaning of Theorem [H 

From above, it can be seen that given the states to be discriminated, the category of the optimum point is not 
determined: it may be a non-singular interior point in the critical feasible region, a point on the boundary of the 
critical feasible region, or even a singular point, depending on the prior probabilities. Taking the states in ([T5]) as an 
example, we numerically obtain the optimum points for three different sets of prior probabilities which result in the 
above three different categories of the optimum points. The results are presented in Table [J and the corresponding 
graphics are Fig. [21 Fig. [3] and Fig. [4] respectively. 



IV. ANALYTICAL RELATION BETWEEN THE OPTIMUM SOLUTION AND THE STATES TO BE 

DISCRIMINATED 

In this section, we shall use the equations established in the previous section to make some efforts on the analytical 
optimum solution to the problem of unambiguously discriminating n pure states. We shall obtain a formula which is 
not a complete analytical solution but can characterize a clear relation between the optimum solution and the states 
to be discriminated, and we shall give an example to show the use of that formula. 



A. Formula 

First of all, we give another form of the matrix X — F in a way similar to Eqs. (2.5) and (2.6) of [Hi. Suppose the 
states l^/^/e) (/c = 1, • • • , n) are represented in an orthonormal basis of the Hilbert space H spanned by {|V^/c)}^^^, then 
each li/jk) has exactly n components in its representation. If the initial state of the system is 1?/^^), then the state after 
the inconclusive measurement result can be chosen as 



|0/c) = VnolV^/c), k = l--'n, 



(50) 




Figure 2: (Color online) The optimum point is a non-singular interior point of the critical feasible region. 




Figure 3: (Color online)The optimum point is on the boundary Br. 



where Ho is defined in Eq. (|6]). Define an n x n matrix C with \(j)k) as its kth column, 

C = (|0i),-- - ,|0„)). 

It can be directly verified that 

X - F = C^C, 

using Eq. ([T]). 

Before presenting the main theorem of this section, we give three lemmas that will be used later as follows. 
Lemma 4- Suppose A is an n x n matrix. Let the adjugate of A be denoted by A*, whose (i^j) entry is the 
cofactor of A. Then 



Rank A* = < 



n if Rank A = n 

1 if Rank A = n-1. 

if Rank A<n-2 
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Figure 4: (Color online) The optimum point is a singular point. 



Lemma 5(Binet-Cauchy). Suppose A and B are matrices of sizes m x n and n x m respectively, then 



det (AB) 





det A-detB 



1 



E detA . 

U<Jl<---<Jrr.<n \jl j 



m 



if m > n 
if m = n 

if m < n 



(54) 



where A 
B 



denotes the m x m submatrix of A whose kth column is the j^th column of A, and 



denotes the m x m submatrix of B whose kth row is the j^th row of B. 

\ 1 • • • ml 

Proof of these two lemmas can be found in many algebra text books. 

Lemma 6. Let C* denote the adjugate of the matrix C defined in Eq. ([5T]) and \ck) denote the transpose of the 
kth. row of C* (/c = 1, • • • , n). If the optimum point p^^^ is a non-singular interior point in the critical feasible region, 
then at p^^^ each \ck) can be written as 



(55) 



where \p) is some normalized vector, ^ is a positive parameter and e*^^ is a phase to be determined. 

Proof. Let Ck denote the submatrix of C by deleting the /cth column \(j)k)^ and {X — T).. denote the submatrix of 
X — F by deleting the zth row and the jth column. Similar to Eq. ([52|) . it can be directly verified that 



so 



(x-rv. =clc„ 



det {X - r), . = det(ClQ) = (-1)*+^' (c.|c,) 



(56) 



(57) 



according to Lemma [5l 

Since p^^^ is a non-singular interior point in the critical feasible region IZcf^ we have Mk (Popt) ¥^ foi" some 
k G {!,••• which implies that Rank (X — F) \p^^^ = n — 1, or Rank C\p^^^ = n — 1, equivalently. Therefore, 
RankC*|p^^^ = 1 according to LemmaSJ implying that all |c/e)'s are proportional to each other. So each \ck) at p^p^ 
can be written as 



(58) 
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where \p) is some normalized vector and ak is a coefficient to be determined. 
From Eq. ([57|) . we have 

Mk (p) = det (X - r)^^ = {ck\ck). (59) 
By substituting Eqs. ([25]) and ([58]) into ([59]), we have 

alak = (60) 

where A > 0. Let A = then Eq. ([55]) holds. 

The factor ^ cannot be zero, otherwise an (p) = would be degenerate at PoptJ contradicting the assumption that 
p^p^ is not a singular point. ■ 

When the states to be discriminated and the prior probabilities are fixed, the difference between any pair of phase 
factors e*^% e^^^ {i 7^ j) is fixed, while the phase factors e*^'^ (/c = 1, • • • , n) themselves can be altered since an arbitrary 
total phase can always be added to all e^^'^'s by choosing an appropriate phase for \p). The phase differences are what 
really matter in the following discussion. 

Theorem 5. If the optimum point p^^^ is a non-singular point in the interior part VLr of the critical feasible region, 
then the components of p^^^ can be written as 

n I 

= e-^^^^e^^^/^(V^,|V^fc), Vz = l,---,n, (61) 
and the optimum average success probability can also be written as 

n 2 



(62) 



Proof. Since p^^^ is a non-singular interior point in the critical feasible region, by substituting Eq. ([55]) in Lemma 
|6]into ([57]), we get 

det(C|Cfc)|p^^^ =(-l)^+'v^e-^'^+^^'^e', Vi,fc=l,--- ,n. (63) 

Noting that the (i, j) entry of X — F is {ijjilijjj) —piSij, the algebraic cofactor of the (i, j) entry is (— 1)*+-^ det (X — F).^. 
and det(X — F) = 0, we perform Laplace expansion on the determinant of X — F along its ith row. 



det(X-r)|p„^^ =^(-ir+'^(V'i|V'fe)det(X-r),fc|p„^^ -PiMi{p,^,) 

n 

= ^(-l)*+'^(V'.|^/'fc)det(ClCfe)|p„^^ -p.M,(p„J (64) 

k=l 
n 

k=l 

for i = 1, • • • , n, where we have used Eqs. ([57|) . ([59]) and (|63|) . Eliminating from both sides of (|64|) . we obtain 

n n 

liPi = Y.(^^\^^^Vin^^~'''^''' = V^^^ Vz,fe = l,... ,n, (65) 

which immediately implies Eq. (|6T]) . Summing up Eq. (|65|) for alH = 1, • • • , n, we eventually get Eq. (|62|) . ■ 

It should be pointed out that a similar result was derived in Ref. |33] for the special case where X is a real matrix, 
and our Theorem [5] can be considered as a generalization of that result to the situation where X is complex. 

It should also be mentioned that Theorem [5] actually gives an analytical relation between the maximum average 
success probability and the n pure states to be discriminated but not a complete analytical solution, since the explicit 
expressions of the phases e*^^ (/c = 1, • • • , n) are not given in Theorem (H However, Theorem [5] may sometimes help to 
simplify the calculation of the optimum solution in special cases, as we shall show in the next subsection. And it may 
also help to obtain some bounds of Popt or work out the phases e*^'' {k = 1, • • • , n) by numerical methods according 
to Remark |3] given later. 
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Remark 2. It seems that the solution (|6T]) and (|62|) would change if the phase of any state \ijjk) is changed. But 
actually the corresponding e*^^ will also be changed in that case and we can see below that any term e^^^ |?/^/e) in Eqs. 
(|6T]) and (|62|) remains unchanged up to a global phase for all e*^^ |?/^/c)'s. 

In fact, if some is transformed as 

\^i)^e'^^i), (66) 

where e*^ is an arbitrary phase while the other n — 1 states stay unchanged, then according to the definitions of \ck) 
and e^^^ in Lemma [H the phase e^^^ changes as follows: 

I e'^^^ (unchanged) if k = i 

Since any global phase can be eliminated from all e^^^'s by the vector Eq. (|67|) is equivalent to 

gi^, ^ fe'^'^ (unchanged) if A; 7^ i 

^ |e^(^fc-x) if k = i' ^ ^ 

Therefore any term e'^^^li/jk) (including the one that k = i) stays unchanged, considering Eq. (|66|) . 

Remark 3. The maximum average probability has an interesting property that it must be the value of a stationary 
point [ill of the expression at the right side of Eq. (|62|) if the phases e*^^ are allowed to change freely. This is because 
Pi is real which requires that 



V7^e"''-(^i|(Ev^e'^N^fe)) - (E^/7^e-''''(V'fe|)^/7^e"'■|^i) =0, (69) 
according to Eq. (|6T1l , and Eq. (|69|l is equivalent to 



d_ " ' 



k=l 



0, (70) 



which is exactly the restriction equation of 6i that must be satisfied when the expression at the right side of 
reaches a stationary point. 



B. Example 



In this subsection, we give an example to show the use of Theorem [5l 

Suppose is a set of linearly independent pure states that (V^ilV^i) 7^ 0, i = 2, • • • , n and {ipiltjjj) = 0, 

Vi,j = 2, • • • , n, let's calculate the maximum average success probability when the optimum point p^p^ is a non- 
singular interior point in the critical feasible region 1Zcf- 

Without loss of generality, we can choose that e*^^ = 1. According to Eq. (|6T]) and the fact that pi is real 
(i = 1, • • • , n), we can directly get 



so 



Substituting Eq. ([711) i^^o 



0k = ATg{{^k\^l))-7T, (/c = 2,- 

Pk = i-^jm^i)\ (fe = 2, 

or using p = ^^^^ ^iPi, we have 



,n) 



Popt 



l-2^V7l7^KV^fc|V^i)l- 



fe=2 



Since < p/c < 1, ^ = 1, • ' ' 7 ^5 according to Eq. ([72]) . we obtain ./^ > |(V^/c|V^i)| (for all /c = 2, 



(71) 



(72) 



(73) 



and 



X]fc=2 Y ^KV^ilV^fc)| ^ 1, which are the conditions for the optimum point p^^^ to be a non-singular interior point in 

the critical feasible region IZcf- 

When n = 2, Eq. (O gives the weh known Ivanovic-Dieks-Peres limit [13, [H, [H, [sO] . 
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V. A GENERALIZED EQUAL-PROBABILITY MEASUREMENT PROBLEM 

A special scheme to discriminate quantum states unambiguously is the so-called equal-probability measurement 
(EPM) [HHl, which requires that the probability of each measurement outcome is equal, i.e., Pi = ■ • ■ = Pn- In this 
section we will treat a generalized version of the EPM problem, which is defined as follows. 

The Generalized EPM Problem (GEPM). If it is required that 

Pi : P2 ' - - ' ' Pn = wi : W2 : - - ' : Wn (74) 

for a given set of non-negative numbers Wi {wi > 0, Vz = 1, • • • , n) when the average success probability of unambigu- 
ously discriminating the states reaches the maximum, one needs to work out the prior probabilities {7i}^^^ 
or the conditions these prior probabilities should satisfy. 

We have the following result on this generalized version of EPM problem using Theorem [3l 

Theorem 6( Generalized EPM). Suppose ^ is a matrix with / ^/wl as its iih. column, let cFmin denote the minimum 
eigenvalue of ^'•'^ and M^^^^ = Mi (p) \pi=wiam,in,,--- ,Pr^^wr^(Trr^ir^ be (n - 1) x (n — 1) principal minor of X — F 
corresponding to its ith diagonal element while pi = wiaminr ' ' iPn = '^nO'min- If M^^^^ > for some i G 
{1, • • • ,n}, then the sufficient and necessary conditions that the GEPM is the optimum POVM to unambiguously 
discriminate the given states is that 

j^GEPM 

li = j^GEPM ' Vz = 1, • • • , n. (75) 

If Mf = for ah i G {1, • • • , n}, then there exist a range of different 7's for which the GEPM is the optimum 
POVM for unambiguous discrimination of the states 

Proof. We suppose that the optimum solution p^^^ is pi = wiT]., • • • ^Pn = WnT] according to Eq. ([74|) . and the 
corresponding F matrix is didig{wir],W2r], • • • ^Wnf]), where rj is to be determined. By some simple calculation, it can 
be shown that in this situation the condition X — F > can be converted to — 7^/ > where / is the identity 
matrix, so (imin, the minimum eigenvalue of is exactly the maximum feasible value of 77, which means that 

Pi — '^i^^min 

(i = 1, • • • , n) is the optimum solution p^^^ when Eq. ([74]) has to be satisfied. If Mf ^ for some 
i G {!,••• 7 ^l? Popt is a singular point, so Eq. ([75]) holds according to Eq. ([25]) . 

On the other hand, if M^^^^ = for alH G {1, • • • , n}, then the point where pi = i^^iCTmm, • • • ,Pn = uJnC^min 
in the critical feasible region is a singular point, so the normal vector changes discontinuously in the neighborhood 
of p = {wiGmin^ ' ' ' ^uJnCTmin) in the Critical feasible region IZcf- Thus, the GEPM is the optimum unambiguous 
discrimination scheme for a range of different 7's. ■ 

The solution of the original EPM problem follows immediately from Theorem [6] by setting wi = - - • = Wn = 1. 

It is obvious that for any set of linearly independent quantum states, it is always possible to find prior probabilities 
{71 7 • • • iln} such that the generalized EPM is the optimum scheme to unambiguously discriminate these states. 

VI. UNAMBIGUOUS DISCRIMINATION OF THREE PURE STATES 

In this section we shall use the results and method presented in Sec. Ill to study the unambiguous discrimination 
problem of three linearly independent pure states, mainly for the non-singular interior optimum solution. 

A. General equations 

Suppose the three states to be discriminated are |V^2), IV^a) with prior probabilities 71,72,73, and they are 
linearly independent. Then according to Eq. ([25]) we can have the following equations 

Ml {p) = (1 -p2)(l -p3) - 10,2^1^ = A71, (76) 
M2 (p) = (1 -pi)(l -P3) - 1(^1 1^3) I' = A72, (77) 



Msip) = (1 -pi)(l -P2) - 1(^1 1^2) I' = A73. (78) 
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From these three equations, we can obtain 



'"^'-V l(^i|^3)P + A72 ' ^''^ 



T „ _ (KV'2|^3)P+A7l)(|(^l|V>3)P + A72) .... 

^'"''i |(Vi|^2)P + A73 • ^''^ 

Substituting Eqs. ([79|) . (|8Q|) and (|8T]) into det (X — F) = and making some rearrangements, we can get the 
fohowing equation of A 



7A^ - SA - 2|Tp + 2v'7A3 + i?A2 + ^A+|T|2Re (T) = 0, (82) 
where Re (T) represents the real part of T and T, 7 are defined as 

7 = 717273, ^ = (V^l|V^2)(^2|^3)(^3|V^l), (83) 



i^ = 7l72|(V^l|V^2)|' +7273|(V^2|V^3)|' +7l73|(V^l|V^3)|', (84) 
S = ll\{M^2)\\^3\^l)?^l2\{^^^^ (85) 

In Ref. it is doubted that whether a general closed form of the maximum average success probability of 

unambiguous state discrimination exists. One can show that Eq. (|82]) can be converted to a polynomial equation 
of degree 6, so it is generally difficult to find an analytical solution for unambiguous discrimination of three states. 
However, in some special situations, analytical solutions can be obtained, and we will give some examples in the next 
two subsections. 

B. Special cases 

For some special cases, Eq. (|82|) can be simplified and one can obtain exact analytical solutions. 
Case 1. Suppose that (V^i|V^2) = 0, but (V^2|'^3) 7^ 0, (V^alV^i) ^ 0, then Eq. (|82|) can be simplified to 

717273^^ - KV^3|V^i)(V^2|V^3)|^73A = 0. (86) 

We can easily obtain that A = or 

^ ^ |(V^3|^l)(^2|^3)| .3^. 
V7172 

and the negative root has been discarded. 

With some observation, it can be verified that A = is not a solution for Eqs. ([79|) - (|8T]) . unless (V^ilV^a) = or 
(V^2|V^3) = 0, which is a trivial case. Thus Eq. (|87|) is the unique solution for this situation. Substituting Eq. (|87|) 
into ([79|) - (|8T]) . we can get the optimum average probability 

p,^, = 1 - 2v^KV^i|V^3)| - 2vW|(V^2|V^3)|. (88) 

This agrees with the result in Sec. lIVBi 

Case 2. Suppose that (V^i|V^2)(V^2|V^3)(V^3|V^i) is purely imaginary, i.e. Re ((V^i|V^2)(V^2|V^3)(V^3|V^i)) = 0, but 
(V^i|^2)(V^2|V^3)(V^3|V^i) ^ 0, then Eq. ^ becomes 

7A^-^A-2|Tp = 0. (89) 
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This is a cubic equation of A, so we can get analytical solutions in general. It can be proved that this equation has 
and only has one positive root of A using Vieta's theorem [69], and this positive root is 

f2772|T|2 + ZJ%l-f*\T\'^ - 37353) ' o 
^ - ^ ^ ' (90) 



37 



where we have used the notations defined in Eqs. (|83|) - (|85|) . 
It can be shown that 



> 277|T|'^ (91) 
by the mean inequality (tO], so Eq. (|9Q|) can be simplified to 

S 

A = 2W — cos-, (92) 

where 



= arccos^yHp. (93) 

Substituting Eq. (|92|) into ([79|) - (|8T]) . we can eventually get pi, p2, Ps and p for the optimum solution. 

Case 3. If ('0i|'02)('02|'03)('03|'0i) is real and non-negative, then it can be directly verified that A = is a solution 
to Eq. (|82|) . which implies that a singular point exists in the critical feasible region. And in this case, Eq. (|82]) can 
be simplified to 

A' (7'A^ - 2-fS\^ - 87|TpA - Q) = 0, (94) 

where in addition to Eqs. (|83|) - (|85|) . another constant Q is defined as 

Q = S^ -m\T\^. (95) 

Since the expression inside the parentheses of Eq. (|94|) is a quadratic polynomial, analytical solutions of A can be 
obtained from (|94|) . But the analytical solution is too complicated to show any practical meaning, so we are not going 
to include it here. 

It should be pointed out that the analytical solutions for three pure states with with real Gram matrix was given 
in Ref. |38j], and our Case 3 generalizes that result. 



C. Equal-Probability Measurement 

In Sec. [Vl we have given the exact solutions of the prior probabilities to the generalized EPM problem. As an 
example, we solve the original EPM problem of three pure states in a direct way. 
Let pi = P2 = Ps = Pepm and substitute it into det (X — F) =0, one can get 



PEPM = 1 ; — 2~3273 ' ^ ' 



(1 - 3i) (^V3y^27(ReT)^- W3 - 9ReT 



2-33 yV2,^27 {KeXf -W^ - 9ReT 
where T is defined in Eq. (j83j| and 



w = \{HiP2)f + lii^imf + liHi^s)^- (97) 

It can be verified that 

27(ReT)^ < 27|T|2 < W3, (98) 
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where the second inequahty can be proved by the mean inequahty ffdl, then Pepm can be reduced to 



3 V3 3 

where 



PEPM = l-2-\l^cos(l-^- ), (99) 



3V3Re T , , 

= arccos — . (100) 

wVw 



Thus according to Eqs. ([76|) - ([78|) . the prior probabihties 71,72,73 must be 

^ |H^COS^(f-f)-|(^2|^3)|^ 

4Wcos2(f-f )-iy 

_ |H^cos^(f-§)-|(^i|^3)|^ . . 

Wcos2(f-f)-Vl/ • ^^^^^ 

^ |M/cOS^(f-f)-|(^l|^2)|^ 

4Wcos2(f -|)-Vl^ 

Remark 4- In Eq. (|10ip . the denominators of 71, 72, 73 can be equal to zero, and in this situation the optimum 
EPM point is a singular point. This can be shown as follows. If the dominators in (|10ip are equal to zero, then 

so 

(9 = or 27r. (103) 

(Note that cos (f — f ) can not be — ^, otherwise Pepm would be larger than 1 according to Eq. (|99|) .) 
Then it can be seen that 

3\/3Re(T) =WVW (104) 

by substituting Eq. (fT03|) into (fTOO|) . 

According to Eq. (|98|) and the condition that " = " holds in a mean inequality (tO], one can have 

\{H^2)\ = \{^im\ = \{H^3)i (105) 

Then according to Eqs. and (|99l) . 

Pi =P2=P3 =PEPM = 1 - KV^l|V^2)|. (106) 

Substituting Eq. (|106p into Eqs. ([76]) -([78]), it can be seen that 

Ml (p) = M2 (p) = Ms (p) = 0. (107) 
So, when the dominators in Eq. (|10ip are equal to zero, the optimum EPM point is a singular point. 



VII. CONCLUSION 



In this paper, we have mainly studied the problem of optimum unambiguous discrimination of n linearly independent 
pure states. We have derived some analytical properties of the optimum solution to this problem, and established 
two sets of new equations in Theorems 3 and 4 which provide detailed methods to obtain the optimum solution 
in different situations. We have also presented a geometrical illustration of the equations we established with a 
numerical example in Sec. IIIIBi An analytical formula which shows the relation between the optimum solution of 
the unambiguous discrimination problem and the n pure states to be identified has been derived in Sec. [iVl And 
we have also solved a generalized EPM problem in Sec. |Vl with the proportion of the occurring probabilities of the 
measurement outcomes to be fixed. Finally, the optimum unambiguous discrimination problem of three pure states 
is studied, and analytical results has been presented for some interesting cases in Sec. [VTl 

It is no doubt that the problem of discriminating quantum states is important in quantum information science since 
it has wide application to quantum cryptography and quantum communication, so it motivates a lot of researchers to 
explore different kinds of optimum discrimination strategies. In addition to many important results mentioned in Sec. 
HI some other novel strategies such as minimax discrimination [65, 66] and maximum confidence discrimination (67l.[68| 
have been introduced recently. We hope that our results presented in this article may stimulate further research to 
the optimum state discrimination problem in general. 
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